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A VECTORIZATI ON OF THE JAMES ON-CAUGHEY NYU 
TRANSONIC SWEPT-WING COMPUTER PROGRAM FL0-22-V1 
FOR THE STAR-IOO COMPUTER 

R. E. Smith, J. I. Pitts, and J. J. Lambiotte 
Langley Research Center 

SUMMARY 

The computer program FLO-22 for analyzing inviscid isentropic transonic 
flow past 3-D swept— wing configurations has been modified to use vector 
operations and run on the STAR— 100 computer. The "vectorized version" described 
herein is called FLO-22— VI. Vector operations have been incorporated into 
Successive Line Over— Relaxation in the transformed horizontal direction. Vector 
relational operations and control vectors are used to implement upwind differ- 
encing at supersonic points. A high speed of computation and extended grid 
domain are characteristics of FLO— 22— VI. The new program is not the optimal 
vectorization of Successive Line Over— Relaxation applied to transonic flow; 
however, it proves that vector operations can readily be implemented to increase 
the computation rate of the algorithm. 

INTRODUCTION 

The program FLO-22 has been developed for analyzing inviscid isentropic 
transonic flow past 3— D swept— wing configurations (ref. 1, 2, and 3)* The 
program originators are Prof. Antony Jameson (New York University) and Prof. 

Dave Caughey (Cornell University) who developed the program under NASA grants 
NGR— 33-016— 167 and NRG— 33-016-201 at the Courant Institute of Mathematical 
Sciencies of NYU. The program is designed to run on third generation computers 
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such as the CDC 6600, and CYBER 175- FLO-22— VI is a first attempt to design 

FLO-22 to run efficiently on the STAR— 100 computer at the Langley Research 
Center. The objectives have been to decrease the computational run time and 
to increase the grid domain. 

The STAR— 100 is a vector processing computer capable of achieving high 
result rates when a high degree of parallelism is present in the computation. 

The STAR— 100 operates on vectors, most simply thought of as arrays in memory, 
as opposed to scalars . When an identical operation is to be performed on 
consecutive elements in memory (i.e. c^=a^*b^ i=l,2, ...N), a vector instruction 
is issued to perform the operations. Each vector instruction involves a time 
penalty, called vector start— up, regardless of the length N of the operation. 

As N increases in size, the operation becomes more efficient since the penalty 
becomes relatively less important. Aronher penalty that generally occurs is 
the necessity for additional storage for temporary vectors instead of temporary 
scalars . 

The goal in designing (or vectorizing a STAR program is to induce long 
vector operations wherever possible. The STAR-100 at the Langley Research 
Center has one half million words of primary memory with virtual memory 
architecture . The operating system automatically brings from secondary memory 
any piece of data which is referred to by the program but does not currently 
reside in primary memory. Hence, overlay or buffering procedures often coded 
into large data base programs can be avoided. More detailed information relative 
to the STAR— 100 computer may be found in references 4 and 5 • 

The algorithm used in FLO-22 is the Successive Line Over— Relaxation technique 
(SLOR) applied with the Murman and Cole differencing scheme (ref. 2) to the full 
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three-dimensional transonic potential equation. There is a coordinate 
transformation (ref. 2) from the physical domain to the computational domain 
which further complicates the nonlinear partial differential equation, tut 
simplifies the handling of the boundary conditions. Basic in the numerical 
procedure of SLOR is the solution of tridiagonal systems of linear algebraic 
equations along lines of the computational grid. In FLO-22 there is the option 
to use lines in the transformed horizontal direction (YSWEEP) or in the 
transformed normal direction (XSWEEP) . Also FLO-22 has the option for grid 
refinement in the solution process. The first step is the computing of a 
solution over a coarse grid. Refined solutions are obtained by calling 
subroutine REFIN which halves the coarse grid and interpolates the coarse 
grid solution for initial guesses for the refined grid. For wing config- 
urations typical of transport aircraft, the SLOR is normally applied in 
the horizontal direction using subroutine YSWEEP with the input variable 
STRIP=1. It is this routine in which the major revisions have been made to 
create FL0-22-V1. The other major revision is the direct addressing of the 
variable G(l,J,K) which contains the potential at all the grid points for 
any given iteration. FLO-22 has been designed to have four spanwise grid 
planes of the potential solution in primary memory at one time. After the 
correction to the potential on the Kth plane has been performed, the Kth+2 
plane is buffered from secondary memory into primary memory and the Kth-2 
plane is buffered out. Since the STAR-100 has virtual memory, FLO-22— VI 
has been designed to address all variables as being in primary memory. This 
has required changes in several subroutines tut primarily MIXFLO and YSWEEP. 
Other changes have been made but do not strongly affect the computational speed . 
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As stated before, the solution of tridiagonal systems of linear algebraic 
equations is the heart of SLOR technique. Typically on the finest grid using 
YSWEEP, over five hundred tridiagonal systems containing 192 linear equations 
are formed and solved per iteration. The efficient solution of a system of 
tridiagonal linear equations is a recursive process and is not in itself 
vectorizable . In addition, since the systems arising from the SLOR method 
are interdependent, they must be solved one at a time. However, in applying 
SLOR to the transformed nonlinear potential equation, the major computational 
activity is the forming of the coefficient matrix and the right hand side of 
the tridiagonal system. These computations are vectorizable, and if the number 
of equations is sufficiently large, the STAR-100 can be efficiently employed. 
Since the largest number of equations are created by sweeping in the horizontal 
direction, the vectorization of YSWEEP is a logical starting point. In FLO— 22— VI 
the coefficient matrix and right hand side of the tridiagonal systems are formed 
along lines in the transformed horizontal direction using vector operations and 
an efficient machine coded scalar tridiagonal solution is performed to obtain 
the correction to the potential. Vector lengths are equal to the number of grid 
points in the horizontal direction. 

A particular problem in vectorizing a program for the solution of transonic 
potential flow is the necessary test for supersonic points and the subsequent 
computation of the upwind differences in the Murman and Cole scheme. Testing 
for supersonic points and modifying the computation inline with a vectorized 
operation for the subsonic points is impossible. 

The problem is overcome at some computational expense by performing the 
supersonic computation with vector operations for all points and by using a 
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control vector to apply the computation only at those points where the flow is 
supersonic . 

The conversion of FLO-22 is complete and several cases have been run on the 
STAR- 100 computer and the solution times compared with those obtained on the 
CYEER 175 • The maximum grid size in FLO-22- VI has been expanded to 500x52x32 as 
compared with 192x24x52 in FLO-22. No serious problems have been encountered. 

This report describes the major routines where changes have been made to 
vectorize FLO-22 to create FLO-22- VI. Comparison of computational times are 
presented. It should be noted that FLO-22— VI is not the optimal vectorization of 
transonic potential flow calculations using SLOR. There is the possibility of 
generating two-dimensional vector operations of parts of the computation in each 
spanwise plane as opposed to each line of the plane. This would give rise to 
longer vectors and higher efficiency. 

VECTORIZATION AND PERFORMANCE EVALUATION OF SLOR IN FLO-22-V1 

The application of SLOR to transonic flow calculations is throughly discussed 
in references 1. 2, and 3> and is not presented herein. Instead, the functional 
operation of the algorithm ir. FLO- 22 ar.d FLO- 22- VI is contrasted. The general 
organization of the program is shown in Figure 1. SLOR is an iterative process, 
and in FLO-22 and FL0-22-V1, each iteration is under the control of the main 
program. After initialization, the main program calls subroutine KIXFLO for each 
iteration. For each plane of grid points in the spanwise direction (Figure 2 ), 
MIXFLO exercises the option for relaxation along horizontal lines using subroutine 
YSWEEP or along vertical lines using subroutine XSWEEP. XSWEEP has not been 
vectorized in FL0-22-V1 and, therefore, will receive no further discussion. Cases 
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requiring XSWEEP will not run with vector operations in FLO-22— VI. Subroutine 
MIXFLO, in addition to exercising the direction of relaxation, performs the 
setup of boundary conditions and in FLO- 22 buffers the grid planes of the 
potential in and out of primary memory. In FL0-22-V1 the potential at any 
grid point is directly addressable. 

Subroutine YSWEEP performs Successive Line Over— Relaxation along each 
horizontal line of the Kth spanwise plane of the grid (see Figure 2), This 
routine has been chosen for vectorization for several reasons. First, FLO-22 
is predominately used at LaRC for transport type wings which require a 
relatively large number of points in the X or horizontal direction, and do not 
require the XSWEEP option . Second, even though the subroutine requires the 
solution to many dependent sets of tridiagonal equations, which is basically 
scalar in nature, it also performs transformation evaluation, PDE coefficient 
evaluation, finite difference approximations and applies the solution of the 
Tridiagonal systems to the potential along each line. Except for the tridiagonal 
solutions, all of the computations are vectorizable . Since the transonic potential 
flow equation is nonlinear and transformed from a physical domain to a computa- 
tional domain, the vectorizable part of the subroutine performs the majority of 
the computation. When the computation for each line on the Kth plane is complete, 
control is returned to l v ! UK FLO. The Murman and Cole difference scheme (ref. 2 and 
3' is implemented in YSWEEP. The scheme changes the difference approximation to 
upwir.d if the flow is supersonic at a grid point. Central difference approxi- 
mations are used for subsonic points. Programmatically, the upwind differencing 
is applied differently in FLO-22 and FL0-22-V1. This will be explained further 


at a later point. 
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In FLO-22, YSWEEP primarily consists of two FORTRAN DO loops. The first 
loop performs the evaluation of all variables through the forward substitution 
in the tridiagonal solution. The second DO loop performs the backward substi- 
tution for the correction to the potential and adds the correction to the 
potential G(l,J,K) along the line. A conditional GO TO sends the control of 
the program to the beginning of the first loop to move through each line of the 
plane. Other peripheral computations are performed but are not discussed here. 
The application of upwind differencing in the Murman and Cole scheme is in the 
form of a comparison of the speed of sound and the latest velocity approximation 
with a branch to a section of the subroutine for upwind differencing if the flow 
is supersonic. This is performed in the first DO loop for each grid point along 
the line. Appendix A shows a listing of the scalar or unvectorized YSWEEP. 

Appendix B shows a listing of the vectorized YSWEEP found in FLO-22— VI. 
Vector FORTRAN instructions and an efficient assembly coded tridiagonal 
equation solver replace the two major DO loops found in FLO-22 and discussed 
above. Using vector operations, arrays replace scalar variables. In the 
unvectorized program all variables for a grid point are computed before going 
to the next grid point. In contrast, a variable in the vectorized program is 
computed for all grid points before going to the next variable. Storage require- 
ments are therefore higher in the vectorized program. In the vectorized YSWEEP 
the computation is not carried through the forward substitution in the solution 
of the tridiagonal equations. Instead, the subdiagonal, superdiagonal, and right 
hand side arrays are computed with vector FORTRAN instructions and stored. This 
is followed by calling a mathematical library subroutine Q4TRIDGE (ref. 5 ) 
for the solution of the tridiagonal system. The solution, which is the correction 
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to the potential at grid points along the line, is added to the potential from 
the previous iteration. In FLO— 22 the test for supersonic points and application 
of upwind differencing is done within the first major DO loop. This computation 
is done for supersonic points only, but this type of testing and branching is 
not feasible for vector computation . In FLO-22— VI a vector relational operator 
establishes a control vector based on the comparison of the latest velocity 
approximation and speed of sound. The control vector is a bit string with a 
"one" stored when the test is true and a "zero" stored when the test is false. 

It is done for all points in FLO-22— VI. The control vector is used to decide 
whether or not to apply the upwind computation in the tridiagonal system. There 
is some loss of efficiency due to the need for extra computations which are not 
used, tut it is considerably faster than reverting to scalar computation. A 
control vector is also used to apply the proper direction of the upwind differ- 
encing along a horizontal line of the computational domain. 

The control in the vectorized YSWEEP is the same as that in the unvectorized 
version. In fact, the vectorization in FLO-22— VI follows the logic of the 
unvectorized program very closely. FLO-22— VI, however, is not the optimal 
vectorization of Successive Line Over— Relaxation for transonic flow calculations. 
As an example, further efficiency can be implemented in YSWEEP by computing 
with vector operations all those variables which are not updated on a grid plane. 
This would mean longer vectors are generated and, therefore, high computational 
efficiency. There is, nevertheless, a trade off between vector length and storage 
requirements . 

Several changes have teen made in other routines in FLO-22 to create 
FLC— 22— ."1 • These changes, for the most part, have minor effect on the 
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computational speed. These changes have been either for primary memory 
addressing of the potential, increased dimensions, or convenience of vector 
operations. Except for incore addressing, subroutine REFIN has required the 
most attention. The interpolation in this routine, combined with buffering of 
planes of the potential in and out of primary memory, has required considerable 
modification; however, vector instructions have simplified the problem (see 
Appendix C ) . 

In FLO-22 the maximum number of grid points for the finest grid solution 
is 192. In FLO-22— VI this has been expanded to 300 . Further, the maximum 
number of points in the transformed wing normal direction or vertical direction 
has been expanded from 24 to 32. This implies that considerably more refined 
solutions can be obtained with FLO— 22— VI. 

The efficiency of FLO-22— VI has been observed by running 3—D wing cases 
on the STAR— 100 computer using FLO-22— VI and comparing the solutions and solution 
times with the same cases run on The CYBER 175 computer using FLO-22. The wing 
configurations that have been used in the case runs are the NACA 65 -jA 012 (stream 
wise) wing at mach .91 and the ARW-2 wing at mach .8. When the same input was 
submitted to both FLO— 22 and FLC— 22— VI, the output results were the same within 
the relative accuracy of the word lengths of the CYBER and STAR computers . Runs 
have also been made for the two wir.gs with the expanded grid in FLO-22— VI. 

Sizing and timing information for the NACA 65^A012 wing at mach .$1 are 
presented in table 1. The table shows the program and computer, the number of 
grid points in each coordinate direction, the total number of grid points, the 
vector length, the central processor time per grid point, and the ratio of the 
central processor time per grid point of FLO-22 /FLO-22— VI . For the vectorized 
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program it is observed that the time per grid point in approximately kO/o less 
when the vector length is 300 as compared to when it is 96 • For the scalar 
program the time per grid point is the same regardless of the number of grid 
points. Another observation from the table is that FLO— 22— VI runs on the 
STAR— 100 1.9 to 3-1 times faster than FLO-22 on the CYBER 175 for the same 
problem solution. With the dimensions of the potential equal to 300x32x32 
FLCW-22— VI requires 464,246 decimal storage locations. 

CONCLUDING REMARKS 

The computer program (FLO-22) for determining 3— D transonic flow past 
swept— wing has been modified to use vector operations and run on the STAR— 100 
computer. The first vectorized version is called FLO-22— VI. Most of the 
vector operations have been incorporated into the subroutine YSWEEP for the 
application of Successive Line Over— Relaxation in the transformed horizontal 
direction. Vector relational operations and control vectors are used to 
implement upwind differencing at supersonic points. Tridiagcral systems of 
linear equations are solved with an efficient scalar assembly coded library 
subroutine. The potential S(l,J,K) is directly addressable as opposed to 
the buffering in and out of primary memory planes of the perennial. 

Increased speed of computation and extended grid domains are character- 
istics of FLO— 22— VI. The increased speed of computation is a function of 
vector length used in the vector operations. In FLC-22— VI me -ector length 
is equal to the number of grid points in the transformed horizontal direction. 
For the cases that have been run on the STAR-100 where the -ector lengths 
have been varied from 96 — 300 , the computational time per grid point has 
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decreased by a factor 1.8— 3-1 as compared to the computational time per grid 
point run on the CYBER 175* 

FLO-22— VI is not the optimal vectorization of Successive Line Over- 
Relaxation applied to transonic flow; however, it demonstrates that vector 
operations can readily be implemented to increase the computation rate of 
the algorithm. 
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I APPENDIX-A 

! LISTING OF SUBROUTINE YSWEEP IN FLO-22 


SUBROUTINE YSWEEP 


COMMON 


COMMON/FLO/ 

COMMON/SWP/ 


G(193,26,<t),S0(193,35),E0(131>, Z0C131), 

IV (193, 35), ITE1(35),ITE2(35), 
AO(193),Al(193),A2(193),A3(193), 

B0(26), 61(26), B2(26),B3(26), 
Z(35),C1(35),C2(35),C3(35), 

XC(35),XZ(35),XZZ(35),YC(35),YZ(35),YZZ(35), 

NX,NY,NZ,KTE1,KTE2,ISYM,KSYM,SCAL,SCALZ, 

YAW, CYAW,SYAW, ALPHA, CA,SA,FMACH,N1,N2,N3, IQ 
STRIP,P1,P2, P3,BETA,FR,IR, JR,KR,DG,IG, JG,*G,NS 
GK1(193,26),GK2(193,26), 

SX(193),SZ(193),SXX(193),SXZ(193),SZZ(193), 
R0(193),R1(193),C (193), 0(193), 

G10(26), 620(26), G30(26), 64 0(26), 61(26), G2(26), 
I1,I2,K,L,N0,LX,MX,KY,MY,T1,AA0,Q1,Q2,TYAW,S1 


J1 -2 

IF (FMACH.GE.l. ) 

C (11-1) « 0. 

DUl-l) ■ 0. 

DO 12 1-11,12 


J1 


R 0 ( I ) 

R 1 1 1 ) 

GK1 ( I, 1 ) 

12 GK1(I,J1-1> 
J 

I 3 
31 BC 


1. 

1. 


« G( I, 1, L) 

- G(I,J1-1,L) 

J1 

12 

-Tl + BKJ ) *C 1 ( K ) 


DD 32 1-11 

AB 

AC 

YP 

A 

H 

FH 

P 

0 

A 

B 

AZ 

BZ 

CZ 

DZ 

DGI 

OGJ 

DGK 


,13 

-T 1 *A 1 (I )*B 1 ( J) 

Tl+Al ( I ) *C 1 (K ) 

SC ( 1, K ) +9 0 ( o ) 

1 * -RQ ( I ) +A 0 ( I ) * AO ( I ) + YP*YP 

R 0 ( L ) /A 
RO ( I ) *A 

A 0 (I)*( 4 »*YP*YP -FH) 

YP*( 4 »*A 0 (I ) + AO ( I ) -FH) 
XZ(K)*XZ(K) -YZ(K)*YZ(K) 

( X Z ( K ) +XZ(K))*YZ(K) 

- AO ( I ) * X Z ( K ) — YP*YZ(lO 

— AO (I ) * Y Z ( K ) +YP*XZ(K) 

H*H* ( P*A -Q*B) -AO ( I ) +X ZZ t K ) 
H*H* ( Q*A + P*B) -AO ( I ) * YZ Z ( K ) 

G ( 1 + 1 , J,L ) - 6 ( 1 - 1 , J,L) 

G ( I, J + 1 , L ) -GM(IiJ-l) 

G ( I , J , L + l ) -GK 1 ( I , J ) 


-YP*Y2ZU ) 
+ YP + XZZ (K ) 
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OGII 

- G ( I +1 > J » L ) 
♦A3 ( I ) *DGI 

-6 ( I » J # L ) 

"G( 1 1 J t L ) 

♦G(I-1, J,LI 

DGJJ 

» G(I#J+1;L) 
-B3( J )*DGJ 

-G ( I / J# L ) 

-G(I* J*L) 

+ G ( I > J-l* L ) 


DGKK. 

1 

DGIJ 

1 

DGiK 

1 

DGJK 

1 

GX 

GY 

U 

V 
W I 

1 

AU 
AV 
GX Y 
GG 
AA 
riZ 
F XX 
F YY 
F XY 
BV 
UU 
VV 

ww 

uv 

uw 

V w 

A X X 
AZZ 
A X Z 
R 

1 

2 

3 

A 

5 


G ( l > J * l + l ) -G( I # J# L ) -G(I,J,L) ♦G(I,J,L-1> 

♦ C3 IK ) *DGK 

G ( 1 + 1 > J + l> l ) -G(I-1>JU>U 

-G< I+l, J-1>L) +G(I-1* J-1#L) 

G( I+l>J>t*l) -6II*iiJ>L-l) 

-6II-1,J,U1) +G(I-1»J#L-1) 

G 1 1* J +1/L+1 ) -GII#J-1#L+1) 

-G ( If J + l * L-l ) ♦G(I#J-1»L-1) 

Al< I)*DGI 
-B 1 ( J )*DGJ 

GX -SX ( I ) *GY +CA*AO(I) +SA*YP 
GY +SA*AG ( I ) -CA*YP 
R 0 ( I ) * ( Cl ( K ) *DGK -SZ(I)*GY +SYAW 

♦ C A*XZ ( K J *$A*YZ(K) + H*( U* AZ +V*BZJ» 

0 + W*AZ 

V *W*BZ 
HMU*U + V*V ) 

QXY + W*W 
D IM ( A A0> • 2*QQ ) 

A Z * S X ( I ) -BZ +FH*SZ(I) 

1. +H*AZ*AZ 

1. ♦SX(I)*$X(I) ♦ H'MZ*HZ 

SX(I) *H* A Z* hZ 

AV — AU* S X ( I ) -FH* W* $ Z ( I ) 

H* AU* AU 
H*8 V*B V 
F H* W* W 
H* AU* BV 
AU*W 
BV*W 

Rl( I)*(FXX*AA -UU) 

FH* AA -WW 

( RO ( I ) ♦PCr<Ij) * ( AZ*AA -UW ) 

— ( AXX*SXX ( I ) ♦AZZ*SZZ(I> +AXZ*SXZ ( I ) ) *GY 

♦ T 1*( AA* ( C Z * GX M DZ-S X ( I ) *C Z ) *GY ) * RO C I ) 

— H* ( C A * ( AU* A U — AV* A V ) MSA +$A)*AU*AV 
-QX Y * ( b * AO ( I ) ♦ V * Y P 

•MW *W ) *( AO( X ) *AZ +YP *B Z ) ) ) 

— WW* ( CA + XZZ ( K ) ♦ SA*Y ZZ ( K ) ) -W*W*(U*CZ +V*DZ)) 
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AXT 

• ABS l AU*AH I ) ) 


AYT 

• ABS ( BV*B1 ( J ) ) 


AZT 

■ AB$(FH*W*C1 IK)) 


A 

- RO( I)*BETA*AA/AMAX1 1 AXT> AYT>AZT, (1 . 

— RO ( I ) ) ) 

AXT 

• A*AXT 

AYT 

■ A*AYT 


AZT 

» A*AZT 


IF (QG.GF.AA) GO TO 33 


AXX 

• AXX+A2 ( I ) 


AY Y 

• ( F YY* AA -VV>* B2(J) 


AZZ 

» AZZ*C2(K> 


AX Y 

• -Rid )*(FXY*AA +UV)*(A8 + AB ) 


A X Z 

• AXZ+AC 


AYZ 

» -RO C I > * (HZ*AA + VW ) *<BC +BC ) 


BP 

■ AXX 


BM 

• AXX 


B 

• -AXX -AXX -Q1MAYY + AZZ ) 


R 

GO TO 35 

■ AXX*DSI I +AY Y*DG J J +AZZ*DGKK 

♦AXY+DGI J ♦AYZ*UGJK ♦AXZ*DGIK + R 


NS 

• NS +1 


S 

- SIGN( I.,U) 


IM 

• I - 1 F I X ( S ) 


IMP. 

■ IM -IFIX(S) 


AXX 

• UU* A2 ( I ) 


A Y Y 

■ VV*B2tJ) 


AZZ 

» RW*C2(K.) 


A X Y 

» 6 • *S*UV* AB 


A X Z 

» 8 • *S + U1n*AC 


AYZ 

■ 8 * *VW ^BC 


6 X X 

• ( FXX*QG -UU ) * A2 ( I ) 


BYY 

» ( f YY + QO -VV)*B 2tJ) 


5 11 

rn ( Fri*QO — WW)*C2lK) 


BXY 

» -( FXY*Uu +UV ) M AB + AB ) 


3 X Z 

« (AZ*QO — U W > ♦ ( AC ♦AC) 


3 YZ 

» -(HZ*QO ♦ ViOMBI +BC ) 


AQ 

« AA/ 30 


DSLTAG 

• 3XX*0GII +BYY + DGJJ +BZZ*DGKK 
♦ B X Y*DGI J +BYZ*DGJK ♦BXZ*DGIK 


DGII 

• G ( I; J * L ) — G ( I M> J > L ) -G ( I M> J > L ) +G( MM* J» L ) 

+ A3 ( I ) *DG I 

DGJ J 

* G(I#J>1) -G(I,J-1<L) -G ( I> J-l» L ) 

-B3 ( J ) *L»G J 

+ GK j. ( I > J -2 ) 

DGKK 

- G(I»J»U -G( I> J# L-l ) -G(I>J»L-1 ) 

♦C3 ( K ) *DGK 

♦ GK2 ( I* J ) 
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DGIJ 

1 

DGIK 

1 

DGJK 

1 

GSS 

1 

8 

BP 

BM 

B 

1 

R 


G(I*J*l) -G(IN,J,L> 

-G ( I * J-l* L ) ♦G(IM*J-1*L> 

G ( I* J* L ) -6 ( I * J * L-l J 

-G(IM*J*L> ♦G(IM*J*L-1) 

G ( I * J > L ) -G 1 1 # J * L-l ) 

-G ( I * J -1 > L ) +G(I,J-1,L-1) 

AXX*DGII +AYY*DG J J *AZZ*DGKK 
+AXY*DGI J +AYZ*DGJK +AXZ*DGIK 

• 5* ( A Q -1 • ) * ( AX X +AXX + AXY +AXZ) 

AQ*BXX -(1. -S)*B 

AQ*BX X -(1. +S ) *8 

-AQ* ( BXX + BX X + Q2M8YY +8ZZ)) 

♦ ( AQ -1. )*(2.*< AXX ♦AY Y +AZZ) + AXY + AY Z 
(AO -l.)*GSS +AQ*QELTAG *R 


35 


37 


IF 

FR 

IR 

JR 

KR 

R 


( ABS (R ) 


1 


32 


A3 


A2 


51 


B 

B M 
3 

C ( I ) 

Dll) 

CG 

I 

DO A 2 M » I 1 
CG 

IF (ABS ( C G ) 
DG 
IG 
JG 
K G 

GK21 l, J ) 

GK1 ( It J ) 

G ( I, J*L ) 

I 
J 

IF 
IF 
IF 


ABS <FR ) ) GO TO 37 


-A YT* ( GK 1 ( I * J -1 ) -G( I* J-l* L J ) 

-AZT*(Gn 1(I»J) -G ( I * J * L-l ) ) 

-A X T -AYT -AZT 
♦ AxT 

-BK*C ( 1-1 ) ) 


LE 
R 
I 
J 
K 
R 

B 

BM 
1./ (B 

b*bp 

8* ( R 
C. 

13 

13 

Dll) 
.LE 

CG 
I 


-BM + DU-1) ) 


-C ( I )*CG 


■ABSIDG) ) GO TO A3 


IJ -KY) 


GKlIhJ) 
G ( l, J , L ) 
G< I* J#L ) 
I -1 
J +1 
31*51*61 


-CG 


II2.GT.ITE21K) ) 13 « ITE2(k> 
( ITE2 (K) .EQ.MX ) 13 * LX 


♦ AXZ ) 



appendix-^ 


52 


61 


DO 52 
LV 

RO ( 1 ) 
Rim 
GO TO 


I»I1* 13 

I ABS ( 1 


31 


62 I-IO» 


62 


N 
I 

IF 
10 
DO 
A 
H 

FH 
AZ 
BZ 
HZ 
F YY 
FXY 
DGI 
DGK 
V 
U 
W 

G(1,KY+1,L) 


I 

1 F 
E 

NO 

EO (NO ) 
N 


( 10. Nt 


71 

IF 

(I.LE.Il) 


I 

9 


E 

9 


IF 

( I V ( I > K ) . 


ZZ 

* 

73 

IF 

(ZZ.GE.ZC 


N 

m 


GO 

TO 73 

75 

R 

9 


E 

9 

77 

M 

9 


G( 1 

,KY+1, L ) 


— I ABS ( I V ( I j K) )) 
AMINO(LV*IABS(IV( I * K ) ) ) 
LV 


♦ SO < I # K ) * SO ( I # K ) 


NO 

LX +1 

(K.LT.KTE1.0R.K.GT.KTE2) GO TO 71 
NX +2 -13 

13 

1. -RO ( I ) ♦ AO ( I ) * AO ( I ) 

RO ( I ) / A 
RO ( I ) *A 

-AO(I>*XZtK> -SO ( I > K ) *YZ (K ) 

-AO ( I ) * Y Z ( K ) ♦SOII#K )*XZ (K ) 

AZASX(I) -BZ +FH*S Z ( I ) 

1. ♦SX(I)*SX(1) +H*HZ*HZ 
SX(I) +H*A Z* HZ 
G ( I + 1 j KY * L ) -G(I-1#KY*L) 

G ( I * KY • L + l ) -GK2<I.KYJ 
SA*AO(I) -CA*SO(I,K) 

Aim*DGI +C A* AO ( I ) +SA*SO(I,K) 

C 1 ( K ) *DGK +SYAW ♦CA*XZ(K) *SA*YZ(K) 

■ G(I#KY-1,L> 

+ < V*( 1 • -H*BZ *hZ ) -U*FXY -W *HZ ) / < F Y Y* B 1 ( KY )) 
IC 

ITE1 (< ) ) GO TO 71 

G(I3#KY#L) -GUO/KY#l) 

NO +1 

EO(NO) +P3*(E -EO(NO)} 

NO 

RETURN 
I -1 

0 » 

NE.l) GO TO 77 

Z ( K ) -TYArfMXC(K) ♦51*A0(I)*AG(1 )) 

( N-l ) ) C- j TO 75 
N -1 


(ZZ -ZO(N-l) )/( ZO(N) -Z 0 { N- 1 ) ) 
R*EO(N) +(1. -R ) * EO ( N-l ) 

NX +2 -I 

-E 
+ E 


G (M,KY+1,L ) 
GK2(rt,KY) 
GK1(H,KY) 
G(H,KY,L> 

GO TO 71 
END 


G(M#KY— 1*L) 
G ( I > K Y-l> L ) 
GK1 ( n , KY ) 

G ( M > K Y > L ) 

G ( 1 > KY> L ) 
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LISTING OF SUBROUTINE YSWEEP IN FLO-22-V1 


SUBROUTINE YSWEEP 


COMMON 

1 

2 

3 

4 

5 
fe 
7 

COMMJN/FLO/ 

COMMuN/SWP/ 

1 

2 

3 


G(301;34;35)>Sd(301;35)#E0(131)>Z0(131)# 
IV(301,35>, ITE1(35)#ITE2(35), 

AO( 301 >* A1 (301)» A2( 301 ) > A3 ( 30 1 ) , 

B0(3A)>B1(3A)#B2«3<»>,B3<34>> 

Z(35),C1(35)>C2(35)>C3<35), 

XC(35)fXZ(35)#XZZ(35)#YC(35)#YZ(35)#Y2Z(35)» 
NX,NY>NZ,KTE1,KTE2, I S YM* KS YM, SC AL > S C AL Z > 
YAW»CYAW>$YAW> A l PHA> C A> S A, F MACH# Nl> N2 > N3, 10 
STRIP#P1#P2# P3#BETA#FR# IR# JR# KR#DG# IG# J G# KG, NS 
GK1(301>3^)>GK2(301*3*»)> 

SX <301)#SZ(301)# SXX (301 )# SXZ( 301)#SZZ<301 )# 

R0< 301)#R1(3C1),C<301)>D(301)# 

G10(3A)# G20(3A)#G30(3A)#GA0(3<»)#&1< 34)# G2<34)# 


4 II# I2#K#L#N0#LX#MX#KY#MY# T1#AA0# Ql# 02, T YAW* SI 

COMMON /VS TAR/ YP(301)#Aj(301)#H( 301 ) > FH { 301 ) > AZ ( 301 ) #BZ(301># 

1CZ <301)#PZ(3G1)#DGI (3Cl)#DGJ<301)#DGK<3 0i)#b(3C'l)#V<3 01)# W(301># 
2AJ<301)#AV(301>#GJ(301)#AA(301)#hZ(301>#FXX(301)#FYY(301),FXY(301) 
3#3V13G1)#AC(3C1 )# JU(3C1)# VV(301 )*Wta(301)f UV ( 301 ) # 0 * ( 301 )#VW(301)# 
4AaX(3C1)#AYY(301)#AZZ1301)#AXZ(3o1)#AXY(301)#AYZ(301)#AXT<301)# 
5AYT(3Ci)#A2T(301 ) # R ( 30 1 ) # DG 1 1 ( 301 ) # DG J J ( 3 01 ) # DGKn ( 30 1 ) # D G I J ( 301 ) # 
6DGIK(301)#DGJK(301)#$UBD(3Q1)#SUPD(301)#DIAG(3G1)#S(301)# 
7DELTAG(301)#H0LD1{301)#HQLD2 (301) #H0LD3( 301) #HGLD4< 301) # HOLDS (301) 
6 > TtMP (301 )>TEM p l(301)»TEMP2(301 ) 

v #$AVE(301)»LU(301)#ITEMP(201)#CG( 30 1) #riOL 0(301) 

BIT d 1 1 ( 301 ) # 3 12 ( 301 ) 

DESCRIPTOR DR 


J1 « 2 

IF (FhACh.oE.l.) J 1 « 3 


C ( ll-l)-C.O 
D ( 11-1) -0.0 
IV- 12 -Il-H 
RO(IliLV) -1.0 
Ri(IijLV)- 1.0 
GKl(ll#ljLV) - G(Il#l#liLV) 

G< 1 1 II# Jl-liLV)- G(Il#Jl-l#LjL\/) 


J- Ji 


13- 12 
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31 LV- 13-11 -*>1 

YPU;LV)»SO(IljKilV> ♦ BO ( J ) 

SAVE(ljLV) -1 *0-R0 ( 1 1 ; L V ) 

TEMP2(1;LV>»YP(1#LV)* YPdjLV) 

TEMPd;LV>- AO( IljLV)* AO(IljLV) 

AJUjLV) » SAVE<ljLV>+ TEMPU)LV) + TEMP2(1)LV) 

H( 1 ; L V > -ROdl»LV)/AJ(l;LV) 

FH(1JLV)«R0( I1,*LV)*AJ (l;iV) 

TEMPI (ljLV)-AOdljLV)* (4.0*TEMP2(ljLV)-FHdJlV) ) 

TEMP2<ljLV) •YP(liLVJ*(4,0*TEMP(l; LV)-FH(1 jlV) I 
A- XZ(K)*XZ(K)-YZ(K)*YZ(K) 

B » ( XZ (K ) +XZ ( K ) J* YZ (K) 

A 2 C 1 J LV ) » -A(MIljLV>*XZ(K) -YP ( 1 J L V > * YZ ( K ) 

B Z ( 1 ; L V ) • -AOdlJ LV)*YZ (K) + YP d » LV )*XZ (K ) 

TEMPUjLV)- H(ljLV>* HtljLV) 

C 2 ( 1 ; L V ) - TEMP(ljLV)*(TEMPl( 1 ) L V ) * A -TEMP21 lj LV ) *B ) 

1 -AOdljLV ) *X Z Z ( K ) -YPdJLV)*YZZ(K) 

DZUjLV)- TEMPUjLV)*(TEMP2djLV)*A + TE MP1 U ) L V ) *B ) 

1 -AO(IljLV)*YZZlK) + YP d J L V ) *X ZZ ( K ) 

DGI lljlV)»6(Ii+b J. L ) L V ) - G d 1-1 » J « L j L V ) 

DGJdJLV)-Gdl> J + lHilV) - GK1 ( II# J-l j L V ) 
DGK(ljLV)-Gdld>L+ljlV) - GK1(I1»J{LV) 

TEMPl(ljLV)- AldljLV) * DGI(lilV) 

TEMP2UjlV)» -B1(J> * DGJdjLV) 

ud JLV )- TEMPI d)LV)-5XdljLV)*TEMP2djLV)+CA*AO(IljLV) 

1 + S A * YP(liLV) 

V(ULV)« TEMP2( ljLV) + SA *AO ( 1 1 J L V )-C A* YPdjLV) 

W d JL V ) ■ ROdljLV)*<Cl<K)*DGKUjlV)-SZdljLV>*TEMP2(ijLV>+SYAw 

1 + CA*XZ<K)+SA*YZ<K)+HdJLV)*(bdJLV)*AZ(l;LV)+VdjLV> 

2 *BZU JLV) ) ) 

AU(ljLV)- bdjLV)+<m;LV)*AZdjLV) 
a v c l ; L v ) • vu;m«-w(i;m*BZ<iiLV) 

TEMPdjLV)-hdJLV) + (UdjLV)*U(l»LV) + VUjLV)*VdjLV)> 

Q a ( i J L V ) « TEMPdiLV) + W ( 1 J L V ) *W d J L V ) 

HCtDldjLV) » .2 * QC(liLV) 

AA ( l; LV ) • VDIM(AAO>HQLDldjLV)jAAdjLV) ) 

HZ ( 1 ; L V ) » AZdJLV )*SXdljLV> - B Z d J L V )+FHdjLV)+SZ( IljLV) 
FXXUjLV)- 1.0+ H(ljLV) *AZdjLv)*AZ(ljLV> 

FYYdjLV)- 1.0+ SX(I1; L V ) *S X ( IljLV ) + H( 1 J L V ) *H Z ( ljLV)*HZ(ljLV) 
FxYdjLV)- SX ( 1 1; LV ) + HlljLV)*AZd;LV)*HZd;LV) 

BVUJlV)- AVUjLV)-AUdjLV)*SXdl;LV)-FHd;LV)*W<ljLV>*SZUljLV> 
UUCljLV)- H( 1; LV ) + AU(1{LV)* AU(ljLV) 

VV(l;LV)» H ( 1 ; L V ) + 3 V { 1 i L V ) * 8V(liLV) 
kw(ljLV)- Fh ( l ; l v ) * W(l;LV)* w 1 1 ; L V ) 

UV(i;LV)» H ( 1; L V ) * AU I 1 ; L V ) * BV(liLV) 
ukii;lv)» AUdiLV)* w(i;lv> 

VWUJLV)- BV(ljLV ) * VdjLV) _ 
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AXX( l}LV)-«l(IliLV)*( FXXC1JLV)*AA( llLVJ-UU(lJLV) ) 
AZZ(ljLV)-FH(l;LV)*AA(liLV)-«W(l;LV) 

axz(1;lv)-c2.o* ro(iulv) )MAZ< ijm*AA(i;m-uim;m > 

HOLOK 1 } L V ) ■-TEMP 2 ( 1 j L V I ♦ ( AXX { 1; L V) + S XX ( 1 1 J L V ) + AZZ ( 1; LV ) 

1 *SZZdl;LV)+AXZ<l;LV)*SXZdliLV>> 

HOLD2ll;LV>-(AAUilV)MCZ<ljLV>*TEMPlU;m + <DZ<l}LV>-SX(Il;m 

1 *CZdiLV> >*TEMP2djLV)) ) * ROdljLV) 

HQLD3 (1 » L V ) • CA*( AU ( 1 i LV ) *AU ( UL V )-AV ( 1 »LV ) *A V ( li L V ) ) ♦ ( S A+S A ) ♦ 

1 AU(liLV )*AV(1JLV) 

HOLD'tlljLV)* TEMP (ljLV )*(Ud»LV )*ACdljLV)+V( 1;LV)*YP(1»LV) ♦2»0* 
I WdiLV)*(AOdltLV)*AZdiLV)+YPdjLVI+BZdiLV)J> 

HOLDS (l; L V I «-WWd» LV )* (CA*XZZ (K )+$A*YZZ (K ) >-tf (1; LV )*»l (li L V ) * 

1 « U ( 1 ; L V **C ZCHLV1+VC15 LV * 

R(liLV>» HOLDl(ljLV) ♦ TlMHOLLiZ d i LV ) -Hd »L V )* ( HOLD3 d» L V ) - 
1 HOLDMliLV) ) + HOLDS (1 J LV ) ) 

AXTUiLV)* AU(1;LV)*A1 ill* LV 1 
AXTdiLV)* VABS(AXT(liLV)iAXT(liLV) ) 

AYT ( 1 ;L V ) ■ BV(ljLV)* 8l(J) 

AYTdiLV)* VABS(AYrCliLV) iAYT(liLV) ) 

AZT(1;LV)- FH ( i J L V )♦ k(liLV) * C1(K) 

AZTUiLV)- VABS (AZT(liLV) ;aZT(1jLv) ) 

HOLDK 1; L V) « AYTdiLV) 

BI1U;LV)» AXT (1 i L V ) . G E . HOL D 1 ( 1 i L V ) 

HuLDl (i;LV>* GbVC TRL <AXT(1;L V), Bill 1;LV ); HOLDK iiLV) ) 

BIlUjLV)- AZT(liLV).GE.SAVtdiLV) 

SAVc(liLV)* QdVCTRLUZT (iiL V } , BI l(li LV ) iSAVE ( IiLV > ) 

31 1 (liL V) «H0LD1 (li LV ) .Gt. SAVE (IiLV) 

SAVcdiLV)- QfaVCTRL(HOLDldiLV)>BIl(l;LV)iSAVE(l;LV)) 
HQLDid;LV)» ROdxi LV )*BETA* AA (li L V ) / SAV E ( 1 i L V ) 

AXT (IiLV)* AXTdj L V ) ♦ HOLDldiLV) 

AYTdiLV)* AYT ( li LV) * nUL D1 ( 1 i L V ) 

AZT(liLV)- AZT(liLV)* HOLDi'liLV) 

SAVe(xiLV)- G(Il,J,LiLV) + G(Il,j,L;LV) 

DGIi (liLV)-G( I1+1,J,L;LV) -SAVE (IiLV) *G (I 1-1, J, L i L V ) +A3 II 1 i L V > 

1 *DGI(liLV) 

DGjj(l;LV)-Gdl,J+l,LiLV)-BAVE(liLV)+G(Il,j-l,LiLV)-B3(J) 

1 *DG J ( 1; L V ) 

DGKKdiLV)»G(Il,J,L+liLV)-SAVE(liLV)+G(ll,J,L-l;LV>+C3(K) 

1 *DGK (IiLV) 

DGIJdiLV)»G(Il+l,J+l,LiLV)-G(Il-i,J+l,LiLV)-G(Il+l,J-l,L;LV) 

1 +G(I1-1> J-l, LiLV) 

DGIK (1; LV) «G(I1+1> Jt L+liLV)-G(Il+l, J, L-i; LV)-Gdl-i, J,L+i;LV ) 

1 +G d 1— 1 f J # L — 1 J L V ) 

DGJK(li LV)»G(I1# J + l>L + iiLV)-G(Il# j-1»L+1;LV)-G(I1i J+ i>L-liLV) 

1 +G( II, J-l, L— 1 i L V ) 
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AXXdiLV)- AXX(liLV) ♦ A2(IliLV> 

AYYdJLVl- (FYY(l»LV)*AAd»LV)-VV(ljLV))+B2(J) 

AZZ(liLV)- AZZ(ULV)* C2(K) 

TEMPlUiLV)- Al<IljLV>*Bll J)*(-T1) 

TEMP2 ( 1 j LV ) ■ Alt Il;LV)*Cl(K)*Ti 

AXYdjLV)»-RldlJLV)*(FXYdjLV)*AAdjLV)*UVdjLV> ) *2 . 0*TE MP 1 <1 J L V ) 

AXZ(liLV)- AXZ(liLV)* TEMP2(ljLV) 

OUST » -T1 ♦ 81 ( J ) * C1(K) 

AYZ ( 1 JL V ) ■ -R0(IljLV)*(HZ(liLV)*AAdiLV) + VWd;LV))*2.0*C0NST 
SUPDdjLV)- AXX(ljLV) 

SUBD dj LV )• AXXdjLV) , „ % , 

DIAGdjLV )«-2.0*AXX(ljLV) -Q1 *< AY Y ( 1; L V ) ♦ A2 Z ( 1 i LV ) ) 

612(1 ;lv;« og<1; lvj.lt. aa(1;lvj 

NNS • Q8SCNT(BI2(ljLV) ) 

TEMP(ljLV)» AXX(1JLV)*DGII(1JLV) ♦ AYYdiLVJADGJJdjLVJ+AZZdjLV)* 

1 DGKK ( li LV ) + AXY (1 J LV ) *DGI J d »L V ) + AYZ d J L V )*DG JKd »LV ) ♦ 

2 AXZ(liLV)*DGIK(l#LV) ♦ R(1»LV) 

R(liLV) * Q8VCTRL(TEMPdjLV),Bl2(ljLV);RdjLV)) 
ail(ULV)* QQ(ljLV) .GE.AA(liLV) 

NS • Q8SCNT(BI1<1;LV>> +NS 
S<ljLV)»VSIGNd.O>UdiLV)jS(ljLV>) 

AXXdjLV) « UU(ULV)* A 2 ( 1 1 J L V ) 

AYY(ljLV) « VV(ljLV)»B2( J) 

AZZdjLV) - WrtdJLVJ* C 2 ( K ) 

AXY(ljLV) • o.*S( l;LV)*UV(liLV)*T£MPl(liLV) 

AXZ(IJLV) » 8.*S(ljLV)+UW(lfLV)*TEMP2(lfLV) 

AYZdiLV) - 6.+i(XjLV)*VWdjLV)+C0NST 
HuLDldJLV)»(FXX{ ljLV)*QQdjLV)-UUd;LV> ) *A2(I1;LV J 
H0LD2(lJLV)»(FYY(l;LV)*&Q(l;LV)-VV(lJLV)) *t2( J) 

HLJLD3UJLV J«(FH(1 JLV)*QG( ljLV)-WWd JLV) )*C2(K ) 

HOLDAdJLV) — <FXYdjLV)*QQdjLV)+UV<ljLV) )* 2.0*TEMPl(ljLV> 
H0LD5dJLV)«(AZdfLV)*QGdJLV)-UWdlLV))*2.C*TEMP2(i;LV> 
TtMPdJLV)- -(HZ( i;LV )*GC( 1JL V)*V*<1 JLV) )* 2.o*CGNST 
TEMPKljLVJ" 1.0 

TfcMPKljLV)- Q8VCTRL(0u(l;LV)>Bll(ljLV)jTENPl(l;LV)) 

AG(lJLV) • AAdJLVJ/TEMPldjLV) 

DELTAGd jLV) »H JlD1( 1 JLV >*DGII d jLV)*H0LD2d JLV)*DGJ JdJ LV) 

1 +HGLD3(ljLV)*DGKK(l;LV)+ HQL 04 ( 1 ; L V ) * DG I J d ) L V ) 

2 +H0LD5(1JLV)*DGIK(1JLV)+ T EM P ( 1 ; L V ) *D G J K ( 1 ; L V ) 
8I2(1;LV)«S(1;LV) .NE.d.O) 

TEMPl(ljLV) « GIIl-1, J,LJLV) 

TEMPillJLV) » Q3VCTRL(odld»J>LiLV)»BI2(l}LV>iTEMPld;LV)) 

T EMP.2 (i J L V ) - G(Ia-2jJ#LJLV) 

TEF P2 (ljL V) - OaVCTRK G(ll+2> J^L JLV bI2(l;L V) ; TEMP2 ll;L VJ ) 

S AVE( 1; LV ) * G(I1#J»LJLV) - TEMPl(ljLV) 

TEMP (1; LV ) » SAVEdiLV) - T EMP1 d } L V ) +TEM P2 ( 1 J L V ) +A3 ( II J L V ) 

1 *DoI ( 1 j L V ) 
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D6I I ( 1; L V ) • Q8VCTRL(TEMP(1;LV)/BI1(1»LV);DGII(1;LV)) 
TEMP1<1;LV>- G(I1-1/ J-1/L;LV) 

TEMPl(liLV)- Q8VCTRL< 6(11+1/ J-l / L i LV > / BI2 <1 i L V > i TEMPI (1 JLV) ) 
TEMP <1» LV) -SAVEdiLVJ+TEMPldiLV)- G( 11/ J-l/ L; L V ) 

DGIJ(liLV) •Q8VCTRL<TEMP<liLV)#BIl(liLV)iDGIJ<l{LV)) 

T EMP 1 ( 1 ; L V ) - G< 1 1-1 / J / L-l i L V) 

TEMPI (1;LV)- Q8VCTRL (G( II +1/ J/L-1;LV>/ 812(1; LV) i TEMPI <1 JLV) ) 
TEMPdiLV)- SAVE( 1;LV)-G( 11/ J»L-liLVI* TEMPI (1 i L V) 

DGIK ( 1* LV ) « Q8VCTRL(TEMPd/LV)/BIl(l;LV);DGIK(l;LV) ) 

DG J J ( 1» L V J * G ( 1 1* J/LiLV)-2.0*G(Il/ J-l/ L J L V ) +GK1 < II / J-2;LV) 

1 -B3UJ *DGJ(liLV) 

DGKK(liLV)- G ( 1 1/ J* Li L V J-2 .0 *G ( I 1/ J / L-l i L V ) +GK2 ( II / J i L V ) 

1 +C 3 ( K ) *0GK { 1 ; L V) 

DGJK(liLV)- G ( 1 1> J/ LiLV J -G C 1 1 > J/L-1JLV) -G ( 1 1/ J -1# L i L V ) 

1 + 6(11/ J-l/L-liLV) 

TEMPdiLV)- AXX(liLV )*DGII ( 1 5 L V ) + AYY ( 1; LV ) ♦DG J J (1 J LV > 

1 +AZZ(liLV)*DGKK(liLV)+AXY(l;LV)*DGlJ( 1 5 L V ) 

2 +ay2(1;lv)*dgjk (i;lv)+axzu;iv>*dgiku;lv> 

TEMPlll;LV)« AQ(1;LV>-1.0 
TEHP2(liLV)- AXY ( 1 ; L V ) +AXZ(liLV) 

SAVE(liLV)« »5*TEMPldiLV)*(2»0*AXX(ljLV)+TEMP2(liLV)) 
HQLD4CULV)- AQ(1 ;LV )*H0LD1(1;LV) 

HQLD(liLV) ■ HOLDMliLV) - ( 1 . 0-S (li LV ) ) *SAVE ( 1 i L V ) 

SUPD(liLV )«Q8VCTRL( HDLDdiLV )/3IldiLV);SUPDdiLV>) 
HJLDdiLV) « HOLDh ( 1 ; L V ) -( 1 . 0 + S ( 1 i L V ) >*S AVE d i LV ) 
$UBD(liLV)»Q8VCTRL(HULD(liLV)/BIi(l;LV>;SUBD(l;LV) ) 
HOLDdiLV) »-AU(l i L V ) * ( 2 . 0*H0 LD1 (1 i L V ) + Q2 ♦ ( HOL D2 (1 i L V ) +HOL D 3 

1 (1;LV>>> ♦ TEMPld;LV)*(2*0+(AXX(ljLV)+AYYd;LV)+AZZ(l;LV)) 

2 +TEMP2 (i; LV > + AYZ ( li LV ) ) 

DIAG(liLV) -Q8VCTRL (HOLD ( 1 i L V ) / B 1 1 ( 1 i L V )iDIAG(l;LV) ) 
HQLD(1;LV)» TEMPl(l;LV)*TtMP(j,;LV) + A<3(l;LV)*DELTAG(l;LV) + 

1 R ( 1 i L V ) 

R(l;iV)-Q6VCTRL(rlDLD(l;LV)/BIld;LV)/R(liLV)) 

ASSIGN DRzR(liLV) 

CALL 06MAX(X'0A'/0/DR/0/IC0uNT/0/CGNST) 

Ip (ABS(CONST).Lfc.ABS(FR) ) GO TO 37 
FR » CONST 
IR - ICGLNT+I1 
JR » J 
KR ■ K 
CONTINUE 

R(liLV)«R (ljLV)-Af f ( 1 ; L V ) * ( GK 1 { 1 1 / J-l i l V ) -G ( 1 1/ J-l/LiLV) ) 

1 -AZT(ljLV) * (GKldl/ J;LV)-G(Il/j/L-liLV)) 

DIAG(1/LV)-UIAG<1 jLV>-AXT(liLV>-AYT(l;LV)-AZTdiLV) 
SUBD(liLV)-SUBD(liLV)+AXT(ii LV) 

I OPT-O 



APPEND IX-B 


CALL QATRIDGEUV*SUBD#DIAG>$UPD#C>D> IOPT#R) 

CG<Il;LV)»R(ljLV) 

ASSIGN DR# CG < II » L V ) 

CALL Q6MAXtX'0A*#G#DR>0# ICOUNT# 0# CONST ) 

IF (ABS(CONST) .Lfc.ABS(DG) ) GO TO A3 
DG “CONST 
IG • ICOUNT+11 
JG • J 
KG » K 

A3 GK2 ( II# J J LV ) • GK1 ( 1 1# J > L V ) 

GK1(I1#J;LV> »G(I1# J#L#LV> 

G(I1#J#LJLV) • G(Il#J#L#LV) -CG(IljLV) 

J“ J + l 

IF (J-KY) 31# 51#&1 
51 IF (I2.GT.ITE2<K> ) I3*ITE2(K) 

IF (ITE2(K).E0.MX) I3-LX 
LV* 13-11 +1 

ITEMP(1#LV)“ VIA6S(IV(Il#K#LV)iITEMP(l;LV)) 

LU(liLV) • 1 -I TE HP ( 1 ; L V ) 

LJ<1#LV>“ VIABS(LUtllLV)jLJtljLV)) 

3IK1JLV) “ LU(liLV) .LT.ITEMP(liLV) 

ITEMP(i;LV)“ Q8VCTRL(LU<l;LV)#Bll(l#LV);ITEMP(ljLV) ) 

50( 1 1 j L V ) • ITEMP (l# LV) 

RKIliLV) “ LU(lJLV) 

GO TO 31 
oi N »N0 

I • LX +1 

IF (K.LT.KTE1.0R.K.GT.KTE2) GO TO 71 
10 -NX ♦ 2- 13 
LV » 13-10 +1 

SAVE(l;LV)«l«d -ROC IO;LV)+AC(IO»LV)*AO(IO;lV)+SC(IO#K;LV) 

1 * SO(IO#K;LV> 

H(liLV) ■ R0( 1C;IV>/SAV£ (liLV) 

F H { 1 ; L V ) • RO ( 10; L V ) ♦SA VE ( 1 i L V ) 

AZ UJLV)— AO(IOJLV)* XZ(K)-SO(IO#K;LV)*YZ(K) 

5 Z ( l;LV)»-AO(10;LV)* YZ(K)+SO<IO#KjLV)*XZ(K) 

HZ ( 1 ; L V ) “ AZllJLV) * $X(IO»LV)-BZ(l#LV)+FH(i}LV)*SZ(IO#LV) 
TEFP(1»LV)“ H(l;LV)* HZ(ljLV) 

FYY(IJLV) » 1.0 +SX(IOjLV)*SX(IO#LV)+MZ(l#LV)*TFhP(l#LV) 

FXY(iiLV) • SX(IO#LV)+ A Z ( 1 J L V ) * TEMP ( 1 ; LV ) 

OGI(liLV) • G(I0+1#KY#L;LV)-G(I0-1#KY#L;LV) 

C3M1JLV) • G(I0#KY,L+1;LV)-GK2<I0#KY;LV) 

vqjlv) « sa*ao(io;lv) -ca* sc(io#kjlv) 

u < 1; L V) “A1 (IO;lV)*DGI(1;LV)+CA*AO(IO;LV)+Sa*SO(IC#n;lV) 

w(1;lv) »ckk )*ogk auv ) +syaw +ca*xz<k> +sa*yz<k) 

G(10#KY+1#L;LV)« G( io# k Y-1#L ; L V ) + ( V ( l;L V ) * ( 1 .G-BZ ( 1 ;l V ) +TEMP( 1; L V ) 
l )-U(1 ;lv)*fxy(i;lvj-w(1;lv)*hz(1;lv) )/(Fyy(1»lv)*bkky) > 
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I - 10 

IF (IO.NE.ITEKK) ) GO TO 71 
E - G<I3,KY,U-G<I0,KY,U 
NO • NO+1 

EO(NO) - EO ( NO ) + P3ME -EC(NO)) 

N ■ NO 

71 IFd.lt. II) RETURN 
I ■ 1-1 
E • 0. 

I F ( IV( I>K).NE.l) GO TO 77 
III • ZIK) - TYAW*(XC<K) + S1*A0 U ) * AO U )> 
73 IF (ZZZ.GE.ZO(N-l) ) GO TO 75 
N • N-l 
GO TJ 73 

75 CONST -(ZZZ -Z0(N-1))/(Z0(N)-Z0(N-1)) 

E • CONST * EO(N) + ( 1 •O-CQNS T )* EO(N-l) 

77 M ■ NX +2: -I 

G(I»KY+1>L)» G ( K Y-l j L ) -E 
G(h,*Y+l>l)» G(I>KY-l,l) + E 
GK£(M/KY) ■ GK11M,KY) 

GK1(H,KY) « G ( M# K Y» L ) 

G( Mf KY* L ) « G ( I > K Y > L ) +t 
GO TJ 71 
END 
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LISTING OF SUBROUTINE REFIN IN FL0-22-V1 


SUBROUTINE REFIN 


COMMON 


1 

2 

3 

4 

5 

6 
7 


G<301,34*35)*S0( 301*35)*E0(131)*ZQ<131)* 

IV( 301* 35)*ITE1< 35)* ITE2< 35)* 

AC( 301>*A1 1301 >*A2( 301>* A3( 301 )* 

B0(34)>B1(34)*B2<34),S3<34), 

Z(35)*C1(35>,C2(35),C3(35), 

XC(35)*X2(35)*XZZ(35)*YC(35)*YZ(35)*YZZ<35)* 
NX*NY*NZ*KTE1*KTE2* ISYM*KSYM, SCAL* SCALZ* 
YAW*CYAW*$YAW* ALP HA* C A* SA*FMACH* Nl* N2»N3* 10 


MX 

9 

NX 

+ 1 


KY 

9 

NY 

+ 1 


MY 

9 

NY 

+ 2 


MZ 

9 

NZ 

+ 3 


MXO 

9 

NX/2 


♦ 1 

MYC 

9 

NY/2 


+ 2 

MZO 

9 

NZ/2 


+ 2 

I V L « 

301 

* MY 




11 


21 


31 


KK»K 

IF (KSYM.EO 
MZO 
K 

KK »K 

DD 53 K ■ 

J 

J J 

I 

II 

G(1I*JJ*K) 

T 

1 

II 
IF 
J 


l I . G T . 0 ) 


,0) GO 
NZ/ 2 
2 


TO 
+ 3 


11 


KK * MZO 
NY/2 +1 
K Y 
MXO 
MX 

• G(1*J*K) 

I “1 

II -2 
Tl 31 

-1 
-2 

TO 21 


42 


54 

52 

53 


GO 
J 

J J * j J 

IF (J.GT.u) GO 
DO 42 J«1*KY*2 
DO 42 I-2*NX*2 
G ( I * J * K ) • «5*(G(I+1»J*K) 

DO 52 I * 1 * M X 
DO 54 J »2 * N Y * 2 
G 4 I » J * K ) » .5+(G(I* J+1*K) 

G ( I * MY* K ) - 0. 

CONTINUE 


♦ G< 1-1* J*K) ) 


♦G(I*j-1*k) ) 
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IF (KSYM.NE.O) 60 TO 111 
MZM • MZO 
MZST • NZ+1 
GO TO 112 

111 MZM • MZO 
MZST ■ MZ 

112 G ( 1* 1> MZST; I VI ) » G( 1> l» MZM; IVL ) 

IF (MZST.EQ.l) GO TO 113 
MZST* MZST - I 

G ( 1> 1> MZS T; I VI ) • .5 * (6<1»1,MZM|IVL) ♦ G( It 1*MZM-1; IVL ) ) 
MZM • MZM -1 
MZST ■ MZST - 1 
GO TO 112 
113 CONTINUE 

TYAW ■ » SYAW/CYAW 

SI • • 5*SC AL 

NO » KTE1 -1 

EC ( NO ) - 0. 

K « 2 

KK*K 

IF (KSYM.NE.O) KK • K+l 
DC 251 K»KK,MZ 


N 

* NO 


I 

* MXO +1 


IF (N. 

L T . K TE 1. OR • < • GT 

• KT E 2 ) GO TO 231 

Ii 

* ITE1CK) 


12 

* I T E 2 < K ) 


Du 212 

1»I1> 12 


DSI 

« SOII + i,M 

“SG(I-1>K> 

DS K 

• SOI I, K+l ) 

“SC( I»R-1) 

S x 

« A1 ( I ) *DS I 


sz 

• C1(K)*JSK 


DGI 

« G( I + 1 > KY > 

K) -b( I-j.>KY>K) 

DGK 

» G(I»KY,K+1)-G(I>KY,K-1) 

R « 

AMIN0(1»IV(I»K)) 

A • 

1.0 -R + A 0 ( I ) * AO ( I ) +S0< I * K 1 * SO < I,K) 

H * 

R/ A 


Frl « 

R* A 


AZ • 

~AO(I)*XZ(KJ - 

SO(l,K)*YZ<fU 

ez « 

-AO(I)*YZ(iO + 

i C ( I , K ) * X Z ( K ) 

HZ • 

AZ *S X -iZ +FH 

*sz 

F YY 

» 1.0 +SX*SX +H*riZ l * rl Z 

F XY 

■ SX ♦ ri*AZ*riZ 


u ■ 

• Al(I)*DGI 

♦CA*AC(I) +SA*SO(I,K) 

W 

• Cl (K) *DGK 

♦ S YAW +C A *XZ(K) +SA*YZ(R) 

V 

• S A*AO ( I ) 

- CA*SO(I >K ) 
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212 6tI,KY*l/K) • G<I,KY-1*K> 

1 MVM1.0 -H+B2*HZ)-U*FXY-W*HZ) /(FYY + BKKY) ) 

NO • NO +1 

EO(NO) « G(I2,KY,K) -G(I1,KY,K) 

N « NO 

I » II 

IF <K.NE.KTE2.0R.YA*.IE,0.) GO TO 231 
221 I -I +1 

M » NX +2 -1 

NO • NO +1 

EO(NO) • G(M,KY,K) -G(I#KY,K) 

IF (I.LT.MXO) GO TO 221 
I ■ 11 

231 I » I -1 

E >0. 

IF (IV(I>K).NE*1) GO TO 237 

ZZ » Z (K ) -TYAw*( XC (K ) +S1* AO ( I ) *A0 ( I ) ) 

233 IF ( ZZ.GE.ZO(N-l) > GO TO 235 
N - N -1 

GO TO 233 

235 R » (ZZ -ZU(N-l) ) / ( 20 1 N ) -ZOlN-l)) 

t - ft * £ 0 ( N ) Ml. -R)*EC(h-l) 

237 M » NX +2 -I 

G ( I > KY+1»K ) ■ G(M,KY-1>K) -E 
G(^/KY + 1»K) • G ( I > K Y — 1 * K ) +E 
IF <IV(I*K).NE.-1) GO TO 241 

GII>KY,K) » .5*G(I>KY,K-l) + .25MG(I»KY,K4-l)+G(r1»KY»K + in 
IF <IV(I>K+1).LT.1) 

1G<I,KY,K) « »5*G(I>KY>K + 1) + «25*(G(IjKY>K-1)+i»(MjKY> K.-1 ) ) 

G(h»KY>M • G l I > K Y > K ) 

G(1,KY-1,K) ■ .5* ( G ( I / K Y > K ) +G ( I >KY-2 > K ) ) 

G(h/KY-1>K) - .5*(G(M>KY#K) + G ( M» K Y-2 t K ) ) 

241 IF (I.GT.2) GO TO 231 
251 CONTINUE 
261 E 0 < NO *1 ) « 0. 

RETURN 

END 



No. Grid Points 

CPU Time Per 
Grid Point 

STAR 

Vector 

Length 

Ratio of CY 175 
CPU Time to STAR 
CPU Time 

X 

Y 

z 

Total 

CY 175 

STAR 

96 

8 

16 

12280 

1 . 4 X 10 -4 

7 . 5 X 10 -5 

96 

1.87 

192 

16 

32 

O 

00 

CT\ 

l.kxlO^ 

5 .X 10 -5 

192 

2.8 

150 

8 

16 

19200 

l.lj-XlO -11 * 

6. 5xio -5 

150 

2.1 

300 

16 

32 

153600 

lAxlO -4 * 

4 .5xlO -5 

300 

3-1 


Total Time = (Time per grid pt.) x (No . of grid pts . ) X (No. of iterations) 


*This case has not been run on the CYBER computer because of storage and 
total time limitations, however, CPU time per grid point is the same as 
that for smaller cases . 








































Figure 1 (General organization of FL0-22-V1) 
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Figure 2 - Computational Domain 
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